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We discuss the quantum-circuit realization of the state of a nucleon in the scope of simple simmetry groups. 
Explicit algorithms are presented for the preparation of the state of a neutron or a proton as resulting from the 
composition of their quark constituents. We estimate the computational resources required for such a simulation 
and design a photonic network for its implementation. Moreover, we highlight that current work on three- 
body interactions in lattices of interacting qubits, combined with the measurement-based paradigm for quantum 
information processing, may also be suitable for the implementation of these nucleonic spin states. 

5j \ PACS numbers: 03.67.Ac, 14.20.Dh, 42.50.Ex,42.50.Dv, 

D ' The last decade has seen a fervent activity in the simulation of complex quantum phenomena through simple and fully control - 

P-h . lable systems. Examples of quantum simulations are now aplenty: the superfluid-Mott insulator quantum phase transition can 

r — ' be simulated using neutral atoms loading optical lattices [£lj] or arrays of coupled cavities with embedded two-level systems J2]. 

fSJ . Molecular energies have been efficiently computed using quantum algorithms and those of H2 reproduced in a photonic quantum 

simulator yfl . Gravitational black holes, Hawking radiation and the Unruh effect have found fertile ground for their simulation in 

Bose-Einstein condensates Jj], Atoms interacting in special ways with light have been used to theoretically simulate important 

models in condensed-matter theory, such as the Lieb-Liniger one [5], while Dirac's predictions on the Zitterbewegung have been 

reproduced in various contexts IsO]. 

The interest in quantum simulation is manyfold. First, important information on the statistical behavior and properties of 

the simulated mechanisms is gathered by dealing with systems that are more easily manipulated and measured. In this sense, 

the development of quantum computing is gaining further significance as a valuable tool that sheds light on difficult problems 

CP studied from a new (information-theoretical) viewpoint. Second, we are going towards the actual realization of the original 

idea of quantum simulations put forward by Feynman JH]. Yet, this very same concept has considerably evolved since then. 

£\J ■ Deutsch has recently pointed out that " [quantum simulation] would be used for smaller things, not things on a larger scale than 

^ ' a molecule. ..Small molecules and interactions within an atom, subtle differences between different isotopes [...]. And of course 

things on an even smaller scale than that" |3]. Buluta and Nori remarked that quan turn simulators would not only provide new 

results (hardly achievable classically), but also allow to test models whose experimental access is either too expensive or beyond 

^ 1 . the reach of current technology, like string theory iflQIl . 

In this paper we move along such lines and propose a sim ple protocol for mimicking the SU(4) quark model of nucleonic 
spin states based on the digital quantum simulator approach iflQIl . We use qubits to encode the information carried by quarks 
and give a quantum-circuit version of the algorithm needed for the achievement of the spin-up protonic and neutronic states 
as described in Refs. lUlll . Our proposal is economic in terms of resources required and very flexible: we design a photonic 
network based on the use of only six modes, which is well within current experimental capabilities, where the spin-up state of 
a proton or a neutron could be simulated. Moreover, we discuss how an observable emulating the intrinsic nucleonic magnetic 
moment can be measured and we highlight the possibility for implementation in an optical lattice-based quantum simulator. It 
should be emphasized that bare nucleonic states (without including gluons), typical of the first quark models, are not the whole 
picture. Quantum chromodynamics (QCD) is clearly the most rigorous theoretical apparatus for the description of nucleons. 
In spite of their inherent limitations, the use of simple quark models has often been very successful in capturing the general 
nucleonic properties. In fact, these models allow for a very reasonable understanding of many aspects of their spectroscopy 
IU2I1 . and a possible theoretical explanation for such agreement is presented in [130. Our goal here is to make a first step towards 
the simulation of subatomic particles, starting with the spin state of neutrons and protons in the scope these simplified models. 

Quark Model for Protons and Neutrons. - In 1961, Ne'eman and, independently, Gell-Mann and Nishijima introduced a 
classification scheme for light hadrons, now known as the Eightfold way, based on the use of their charge Q and strangeness 
S t . Charge and strangeness are examples of quantum numbers associated to internal symmetries, i.e. symmetries following 
from transformations that do not involve space and time. Other examples are the baryon number B and the isospin /. All such 
quantum numbers are necessary for the understanding of subatomic reactions driven by strong and weak forces. Not all the 
internal quantum numbers used to describe baryons in the Eightfold way are independent for they are bound to the empirical 
Gell-Mann-Nishijima equation Q — h + j(B + S,), where ^3 is the projection of the isospin on an arbitrary axis. 

In line with the Eightfold way, hadrons can be divided in two groups, according to their spin: baryons (fermions) and mesons 
(bosons). Gell-Mann-Nishijima formula can be applied equally to baryons and mesons, thus suggesting the existence of a 
fundamental explanation of the symmetries in the Eightfold way and the range of values that each quantum number can take. 
Such features can indeed be explained by assuming that hadrons are formed by the combination of a small number of constituents 



called quarks. In the first quark model IU4I1 . three types of quarks were proposed to accommodate the known baryons and mesons. 
They are called flavors and commonly denoted as u {up), d {down) and s {strange). They differ from each other in terms of mass 
and quantum numbers. This aspect is illustrated in TableU In this picture, baryons (mesons) are seen as the composite states of 
three quarks (one quark and one antiquark). Modern quark models now include more quark flavors forming heavier hadrons. 

As a first step in the simulation of subatomic particles, here we are interested in nucleons (protons and neutrons), which are 
the lightest particles in the spectrum of baryons. The proton is a uud bound quark state belonging to the symmetric spin-| 
octet (part of the Eightfold way). The neutron also belongs to this octet, but its quark content is udd. From Table [I] one can 
easily check that these quark contents lead directly to the fact that nucleons are indeed baryons with baryonic quantum number 
B=| + | + | = l and that the proton is charged with <2 = | + I — 4 = 1, while the neutron is neutral. The derivation of the 
flavor-spin bound quark states involves the standard addition of angular momenta, as described in detail in Ref. fllll . One is 
eventually led to the following states describing a spin-up proton or neutron 



Iproton) = {\ Ps )\Xs) + \Pa)Wa))I V5, 
Ineutron) = {\n s )\xs) + \n A )\XA))l V2, 



(1) 



where we have introduced the flavor states \p A ) = {\udu)-\duu}) / y2,\ps) - {\udu)+\duu) - 2\uud})/ V6, \n A ) = 
{\dud)-\udd)) I v2, \ns) - {\dud)+\udd) - 2\ddu))/ "Vo and the spin states 

\xa)=Q TIT) - I ITT»/ V2, \xs)=Q UT) + I ITT) - 2| TU))/ VS. (2) 

The subscripts S and A in the states above draw the attention to their symmetry and antisymmetry properties with respect to 
interchange of the first two quarks. The states in Eq. (f2]i are evidently non-separable, thus showing the entanglement in the 
nucleonic state. However, it is important to remark that such entanglement can not be used directly as a resource given the 
impossibility to coherently manipulate individual quarks in the nucleons. For symmetric ground states of baryons, a problem 
arises due to their fermionic nature which enforces their description by overall antisymmetric states. This led to the discovery of 
an additional degree of freedom called color (red, green, or blue), and to the postulate that all hadrons are colorless, i.e. the color 
degree of freedom is described by an antisymmetric state (color singlet). The color structure of all baryons is thus identical and 
so we omit it in the following. However, it is important to emphasize that quarks interact strongly by exchanging color. 
Quantum Circuits.- Let us employ the logical encoding |u)=|0), \d)=\\), |T)=|0) and |J,)=|1). With this notation 

\ PA ) = (eLi&xj) \n A ) = \xa) = (1010) - |100»/ V5, 

(3) 
\ps) = (®j =1 <W)|ns> = \Xs) = (1010) + |100) -2|001»/ Vrj. 

The neutron state can be easily obtained from the equations above simply by applying cr x operators to the qubits embodying 
the flavor degrees of freedom. A detailed construction of the proton state from a fiducial initial state is described in the following 
paragraphs. In what follows, CGy^..)* is used to indicate a generic many-qubit controlled gate G where {jxji..) are the control 
qubits determining the action of the gate on the target k. The notation j l is used when 10)/, activates a controlled gate, while & a j< 
indicates the a-Pauli matrix {a-x, y, z) for qubit k. For easiness of discussion, it is convenient to refer to the constituents of (Q]i 
as qubits j=\, ..,6. Without affecting the generality of our study, we consider the fiducial state \\jjq) - |000) i23|000)456 (other 
choices would result in minor modifications to what follows). It should be noticed that such a choice is unrelated to the physical 
quark model for nucleons, but is just a convenient choice for our protocol. We now apply a Hadamard gate H2=(<t Zj 2+o"x,2)/ v2 
to qubit 2, followed by a controlled-NOT gate CNOT(2)5=|0)2(0|®l5+|l)2(l|®o",t l 5 and look for the unitary transformation < U 
such that 

\ PA ) = #|000>123, \PS)= #1010)123, 

\XA) = #1000)456, \XS) = #1010)456- 

In our scheme, the first three (last three) qubits are clearly used to embody flavor (spin). It is worthwhile to notice that the only 
entangling operation required between flavor and spin qubits is the initial CNOTpjs gate. There is no need to synthesize a full 



TABLE I: u, d and s label the up, down and strange quark. B is the baryon number, Q is the charge, / 3 stands for the projection of isospin on 
an arbitrary axis, while S t is the strangeness. 
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FIG. 1: (a) Circuit for the preparation of Iproton). We show the symbol for H 2 and CNOT (2 )5. The unitary transformation *!/ is decomposed as 
in (b). We show the symbols for single-qubit phase-shift Z=& : , controlled-rotations by angles &, <p and a CCNOT. Here, i?=arccos(- V2/3) 
and <p=n/4. Empty (filled) dots indicate control operated by state |0> (|1)) of the corresponding qubit. The three X = & x gates in the dashed 
box of panel (a) are required for the generation of a neutron state [see Eq. (Hi] . 



six-qubit unitary operation, which represents a remarkable simplification in design due, in part, to our choice for information 
encoding. When the class of allowed gates is restricted to single- and two-qubit ones, any three-qubit unitary operation can be 
implemented with at most 20 CNOT gates and arbitrary single-qubit rotations IU6I1 . In what follows, we provide an explicit quan- 
tum circuit that implements the transformations in Eq. (0| using only 6 single and two-qubit gates. A sketch of the whole protocol 
and the decomposition of < U are given in Fig. [TJ Our procedure is as follows: we concentrate on the triplet (1,2,3), although 
our analysis is applied to (4,5,6) with no changes. The eigenvectors of the reduced density matrix pi23=Tr456(|proton)(proton|) 
all have real components. Therefore, the matrix P that diagonalizes P123 is special orthogonal and unitary and can embody f U. 
However, as the only transformations required by our protocol are those listed in Eq. (0J, there is no necessity for reproducing the 
whole P: as just |000)i23 and |010)i23 are involved, only the first and third columns of < U have to be identical to the homonymous 
in P. This makes the remaining states of the three-qubit computational basis don 't-care logical entries such that { < li-P)\q)m=0 
[^0] for q-Q, 2 [q + 0, 2] (here q-Q, .., 8 corresponds to the binary number identifying the elements of the computational basis). 
We thus see that 



^CR ( 3)2(^)CCNOT 3 5 ZjCNOTj H 2 CR ( 2)3(t?) 



(5) 



where CRy)*(£) is a gate rotating the target qubit k by an angle £ when the control qubit / is in \\)j. Here, (■&, <p) = 
(arccos(- V2/3), n/4). The =* sign in Eq. © is due to the fact that the quantum circuit in Fig. [Tib) transforms |010)i23 (|010>456) 
into —\ps) (— Ifs))- However, this is immaterial to us as \ps)\xs) appears in the protonic state. We stress the key role played by 
CR(2)3(#): for state |000)i23 (|000)456), this is an inactive gate so that the qubit 3 (6) would never be in 1 1 >3 (|1)6). This makes the 
sequence CRp^i/^CCNOT.™ redundant for this entry and < U becomes a Bell-entangling circuit IU5I1 . 

Let us now briefly comment on the resources needed for the simulation of < U. As the rotations used to build up 
Cf\j)k(0 read ^a (^)= sin £ cr^ + cos ^cr^ (£=§,<p), they can be simulated by means of a single CNOT gate and two rota- 
tions each as IU7I1 CR; ; ')a(£)=[1 j ® ^t(£/2)]CNOT(j)jt[ll ; - ® Rk(£/2)]. Moreover, an economic simulation of the Toffoli gate 



CCNOT^j required in our scheme is also possible. Following Ref. 11711 . we can simulate a gate congruent to CCNOT,™ as 

X 2 Wi(f)CNOT (2 )ilVi(f )CNOT (3 )ilV 1 , (f )CNOT ( 2)iW 1 l (f)X 2 , where W k (0 = RkiO^k- This gives a truth table identical to the 
one for CCNOT™! but for the entry corresponding to state 1 1 1 1 }, which are transformed into -| 1 1 1 }. However, as this state does 

not enter into the parts of 'U needed in Eq. dfQ, the congruent-gate simulation fulfills our needs. Our decomposition of < U thus 
requires only 6 CNOT's, far less than the estimated upper bound given in Ref. 11611 . We believe that, although we cannot claim 
for optimality and further improvements may be in order, the decomposition we propose could well be seen as rather efficient. 
The realization of the neutron state goes along the same lines, although it requires three X gates on 1 , 2 and 3 [cfr. Fig. [TJ(a)] . 

The resource-estimate given above is based on the assumption that only two-qubit interactions would be available in the spe- 
cific set-up used for the simulation here at hand. However, in many cases this might well be too limiting as genuine multi-qubit 
interactions could be in order. In NMR, for example, three-body interactions can be efficiently achieved and have been used to 
investigate ground-state properties at criticality 1181 . Analogously, theoretical schemes have been put forward for the achieve- 
ment of genuine and tunable three -body couplings in triangular-cell optical lattices loaded with two-species cold atoms 11911 . 
This scenario is particularly interesting as it offers the possibility for the combination of laser-induced three-body interactions 
and measurement-based quantum computing [20] in a lattice of many-body systems. In this context, indeed, a bowtie lattice-cell 
configuration has been shown to naturally entail three-body couplings suitable for the economic simulation, via a measurement- 
based approach, of CCNOT gates (along with the standard toolbox of CNOT and single-qubit rotations) 12111 . It is interesting to 
notice that both the methods at the basis of Refs. lfl9l 12111 rely on the induction of the couplings "H c - £,- & x j-\& z j& x j+\, among 
others, where i is the label for physical qubits in the neighborhood of a given lattice-cell configuration. This is the so-called clus- 
ter Hamiltonian IU9I1 . whose ground state encodes a cluster state, which is the key resource for measurement-based processing 
of information. Therefore, a very promising setting for the simulation of the nucleonic spin states involves the natural encoding 
of cluster states in lattices of quantum many-body systems coupled by suitable two- or three-particle interactions, as discussed 
in Il22ll . having a lattice-cell structure suitable for the simulation of Toffoli gates. This would allow for the implementation of 




FIG. 2: (a) Protocol for the simulation of a protonic state in an optical network. The single-photon entangled state \4i+)n is generated off-line 
and entangled to modes 5 and 6, both prepared in their vacuum state. Transformation *V completes the protocol. It is realized using three beam 
splitters and a phase shifter, as shown in panel (b) for the mode-triplet (;', j, k) with i=l,4, j=2, 5 and k=3, 6. 



our proposal with just 13 two-qubit gates or, significantly, only 9 two- and three-qubit gates. 

The discussion above does not exhaust all the possibilities for the realization of the nucleonic spin states. In fact, one can 
easily work out a different version of our protocol that is very suitable for an all-optical implementation. Our starting point is 
the fact that Eqs. involve only three elements of the computational basis. The unitary matrix 'ZY (an 8x8 matrix) in Eq. (HJi 
can thus be replaced by the more compact unitary (in the ordered basis {|001}, |010), |100}}) 



<V = 





1/V2 

-1/V2 



-V2/3 1/V? 
1/V6 1/V3 

l/VS i/V3, 



(6) 



where the first two columns are the components of states |?xs) {r—p,x) and the last one is determined by imposing unitarity. 
Eq. © helps in constructing Iproton) if, instead of state \ij/i), we use |*^2)=|00>i4(|0101) + 1 1010))2356/ v2, which can be obtained 
using |</r + >=(|01> + |10»2 3 / V2 as shown in Fig. Ha). 

By encoding the states of the basis used to express *V into spatial modes of light 12311 . we deal with a three-mode unitary 
operation that can be decomposed in terms of beam splitters, phase shifters and rotators as proven in Ref. 112411 . By writ- 
ing the beam splitter operations on spatial modes m and n as BS„,„(w)= cos u>fr x + sin <±>& z , we find < Vy* = SjTtTjT,- with 

T,=BS#(- arccos(l/ V3)), f ; sBS«(-3tt/4), f^=B"S, 7 (0) and Sj being a ;r-phase shift (?'=1,4, j=2, 5 and k=3, 6). Therefore, the 
simple interferometric setting shown in Fig.|2jb) realizes the desired transformation. On the other hand, state |^) should be fully 
within the grasp of current experim ental abilities in bulk-optics and can be implemented starting from a Bell state, as generated 
from a Type-II parametric down-conversion process, and using the handiness of polarization beam splitters as polarization-path 
CNOT gates 02511 . Interesting multi-photon states (of up to six photons) are routinely prepared and manipulated in many op- 
tical labs 12611 . putting our proposal fully within the realm of realistic and implementable schemes for photonic quantum state 
engineering. 

As an application of the emulated quark states described above, we now discuss the measurement of the magnetic mo- 
ments. The net magnetic moment of a nucleon is simply the sum of the moments of the three constituent quarks, i.e. 
//yv=2, = i/ u /(A/ r |o" z> ,|AT) with ,/V={proton, neutron) and //, the magnetic moment of the quarks. By assuming equal mass for 
the d and u quarks, one has jj u - -2//^ ifllfl leading to /Uneuaon //"proton = -2/3, which is in excellent agreement with the value 
determined experimentally. Using our encoding, the nucleon magnetic moment operator reads E- 2^ =1 (|1),(1|-2|0),(0|)<8> cKi+3, 
where the summation is over the constituent quarks. As a only involves local operators of qubits i and i + 3, the measurement of 
(S) will be possible in both the simulation contexts discussed here. 

Conclusions.- We have proposed a quantum circuit for mimicking the nucleonic spin states resulting from the combinations 
of their quark components, as predicted by the Eightfold way. It is important to emphasize that dynamics is not accessible in 
our protocol, which mainly embodies a kinematic model. However, our proposal represents a first attempt to bring the realm 
of quantum simulation to the elementary-particle domain, which was so far unexplored, to the best of our knowledge. Clearly, 
the formulation of a quantum mechanical framework for the simulation of particle physics where also the complex aspects of 
QCD were included would be a very important achievement due to the intrinsic difficulty to simulate it in a classical computer. 
This is an important and clearly difficult topic to be investigated from now on. Here, we moved on in a different direction, and 
presented the simulation of nucleonic spin states as given by a simplified kinematic approach of quarks. We have presented a 
simple generation scheme involving manipulations which are currently implemented in different experimental setups. The most 
demanding part of quantum circuit decomposition proposed here {'U implementation) requires only 6 CNOT's, far less than the 
estimated upper bound given in Ref. IU6I1 for a general three-qubit gate, and we believe that our proposal of decomposition could 
well be seen as rather efficient, although a general proof is still lacking. 

Our scheme does not include color-exchange through gluon emission and absorption and still has room for further develop- 
ments, including future circuit decompositions and physical implementations of other hadron simulations. A possibility is the 
inclusion of heavier quarks or the simulation of mesons. With the addition of strange quarks, for instance, a new scenario would 



be opened, where baryons such as £ + ~'° could be simulated. As a visionary goal, it would be very interesting to design a unified 
simulation framework able to encompass more quark bound states. 
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